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Abstract 

The motor part of a crawling eukaryotic cell can be represented schematically as an active continuum layer. The main 
active processes in this layer are protrusion, originating from non-equilibrium polymerization of actin fibers, con¬ 
traction, induced by myosin molecular motors and attachment due to active bonding of trans-membrane proteins to a 
substrate. All three active mechanisms are regulated by complex signaling pathways involving chemical and mechani¬ 
cal feedback loops whose microscopic functioning is still poorly understood. In this situation, it is instructive to take a 
reverse engineering approach and study a problem of finding the spatial organization of standard active elements inside 
a crawling layer ensuring an optimal cost-performance trade-off. In this paper we assume that (in the range of interest) 
the energetic cost of self-propulsion is velocity independent and adopt, as an optimality criterion, the maximization of 
the overall velocity. We then choose a prototypical setting, formulate the corresponding variational problem and obtain 
a set of bounds suggesting that radically different spatial distributions of adhesive complexes would be optimal de¬ 
pending on the domineering active mechanism of self-propulsion. Thus, for contraction-dominated motility, adhesion 
has to cooperate with ’pullers’ which localize at the trailing edge of the cell, while for protrusion-dominated motility 
it must conspire with ’pushers’ concentrating at the leading edge of the cell. Both types of crawling mechanisms were 
observed experimentally. 


1. Introduction 

Eukariotic cells are spatially extended active bodies that can steadily self-propel in viscous environments at low 
Reynolds numbers [1, 2], It has been understood [3, 4, 5] that in these conditions a combination of stationarity 
and linearity of friction leads to kinematic reversibility and that a symmetric (under time reversal) stroke cannot 
produce self-propulsion. For Stokes swimmers, a variety of non-symmetric motility strategies have been proposed 
and optimized using various efficiency criteria [6, 7, 8, 9, 10, 11, 12]; similar models for crawlers advancing on a 
frictional background were considered in [13, 14, 15]. Most of the self-propulsion mechanisms proposed in these 
papers are fully kinematic in the sense that the time dependence of the shape of a swimmer/crawler is prescribed. This 
implies that appropriately chosen actuators can always perform the required internal movements. In cells the role of 
such actuators is played by active agents and in this paper we focus on the fact that their dynamics, while being driven 
biochemically, must be compatible with the fundamental balances of mass and momentum. 

In the context of eukaryotic cell motility, a prototypical scheme of self-propulsion includes protrusion through 
polymerization of actin filaments which is accompanied by dynamic assembly of focal adhesions; myosin-driven 
contraction of the actin network which allows the motor part to advance a cargo, and, finally, detachment of adhesive 
contacts with a simultaneous depolymerization of actin fibers [16]. It is usually assumed that active polymerization 
ensuring protrusion can be described as the work of spatially distributed pushers, generating positive force couples, 
while active contraction can be viewed as an outcome of the mechanical action of distributed pullers, responsible 
for negative force couples [17, 18, 19, 20, 21]. The role of ATP in reversible adhesion of adhesive patches (focal 
adhesions) is understood rather poorly and they are usually treated as passive viscous binders whose spatial distribution 
may be regulated actively [22] . 
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Since our knowledge of the mechanism controlling the transport and the intensity of active agents performing 
protrusion, contraction and adhesion is rather limited, we adopt in this paper a semi-kinematic approach and treat the 
corresponding distributions as functional control parameters constrained by the fundamental mechanical balances. We 
then pose a variational problem of finding the optimal temporal and spatial distributions of these parameters inside a 
crawling continuum body. In view of some successful attempts to justify such reverse engineering approach [23], we 
anticipate that our optimal solutions will be eventually backed by an appropriate constitutive theory. 

While the real organisms are expected to optimize some measure of a trade-off between the velocity of self¬ 
propulsion and the corresponding energy expenditure, in this paper we make a simplifying assumption that the en¬ 
ergetic cost of self-propulsion is fixed and use as optimality criterion the maximization of the overall velocity. We 
are interested in steady translocation and assume that the internal distributions of mechanical parameters are com¬ 
patible with the traveling wave ansatz. This simplifying assumption allows us to replace the optimization of the 
crawling stroke in space and time by a purely spatial optimization of the internal distribution of active elements in 
the co-moving coordinate system. In the interest of analytic transparency we use the simplest ID model of actively 
contracting continuum (subjected to viscous frictional forces) that has been used repeatedly in the cell motility studies 
[24, 25, 26, 27, 28, 29, 30, 16, 31]. In a setting similar to ours the dependence of cell velocity on the distribution of 
active stresses and adhesion properties was studied in [18] where both contraction and protrusion were represented by 
active couples. Given that protrusion is usually localized at the leading edge of the cell, we model the effect of active 
polymerization differently by using Stefan type boundary conditions on the edges of a crawling segment that fix the 
influx and the outflow of actin, see also [32, 33, 21], For the given strength of protrusion, we prescribe the average 
level of contractile activity, and then search for the optimal internal distribution of contractile and adhesive units. 

Our analysis of the ensuing variational problem demonstrates that radically different distributions of focal ad¬ 
hesions are optimal depending on the domineering active mechanism of self-propulsion. Thus, for contraction- 
dominated motility, focal adhesions have to cooperate with pullers which end up localizing at the trailing edge of 
the cell while for protrusion-dominated motility they must conspire with pushers which concentrate in our model at 
the leading edge of the cell. Both types of crawling mechanisms have been observed experimentally. 

The paper is organized as follows. In Section 2 we formulate the model. A simple analytically tractable case of 
contraction dominated motility is treated in Section 3. In Section 4 we obtain analytically upper and low bounds for 
the self-propulsion velocity in the general case. The same optimization problem was studied numerically in Section 
5. In Section 6 we present some evidence that our lower bound may in fact coincide with the optimal solution. The 
last Section 7 contains the discussion of our results. 

2. The model 

Following [29, 32], we consider a one-dimensional segment of viscous active gel representing the cell lamel- 
lipodium on a frictional substrate. The segment has two free boundaries which we identify as the trailing edge L (t) 
and the leading edge l+(t). The force balance can be written in the form 

d x a = %v (1) 

where cr(x, t) is the stress and v(x, t ) is the velocity. We assume that the frictional coefficient mimicking the distribution 
of focal adhesions [33, 34, 32, 35, 36, 37] is space and time dependent £(x, t) > 0. The constitutive behavior of the 
gel is modeled by the equation [29, 32] 

a = r]d x v +x, (2) 

where the active pre-stress ^(x, f) > 0 accounting for the presence of myosin molecular motors [38, 39, 40, 41, 42] 
is also assumed to be a function of space and time. For simplicity we assume that the bulk viscosity coefficient 
t] > 0 is constant. The assumption of infinite compressibility in (2) allows us to decouple the transport of (actin) 
density p(x, t ) from force balance making the mechanical problem ’statically determinate’. The mass balance equation 
d t p + d x (pv) = 0 can then be solved independently after the velocity field is determined [43, 44]. 

We further assume that some internal mechanism (stiffness of the cell cortex [45, 46, 47, 48, 49, 50], osmotic 
pressure actively controlled by the channels and pumps on the cell membrane [51, 52], etc.) maintains a given size 
Lq = 1+ - L of the cell. Therefore the stress at the edges must be the same cr{L{t), t) = cr(/ + (f), t ) = erg, where o"o(f) 
is then an unknown function. To model active protrusion we impose two kinematic Stefan type boundary conditions 
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characterizing the rate of actin polymerization v+ > 0 and depolymerization v_ 
segment [29, 34, 32, 33, 43] 

'l± = v± + v(l ± (t)J). 


> 0 on the boundaries of the moving 

(3) 


For consistency, the overall mass balance must be also respected on the moving boundaries and we set p(L(t), t)v = 
p(l+(t),t)v+, which implies an instantaneous recycling of depolymerized actin from the trailing edge to the leading 
edge, see [21, 43] for more detail on such closure of the treadmilling cycle. While there is considerable experimental 
evidence that active polymerization is indeed localized at the leading edge of a crawling cell, the de-polymerization 
may be spread along the length of the lamellipodium [33, 32], However, in the interest of analytic transparency, such 
spreading will be ignored in this study, see though [43] . 

The two functions x and £ can be interpreted as infinite dimensional controls parameters and found through an 
optimization procedure. Even in the absence of a detailed microscopic model governing the rearrangement of these 
agents we still need to subject them to integral constraints prescribing the average number of adhesion complexes [42] 


I W+(t) 
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where > 0 is a given constant and 


g(x, t)dx = 


i pLU) 

— I x(x,f)dx=x* 

To Jl_t 


(4) 


(5) 


To Ji_( t ) 

where x* > 0 is another given constant representing the average number of contractile motors [53] . It is clear from 
(4,5) that since we prescribe the density of active agents, the performance of the self-propulsion machinery will be 
proportional to the length of the active segment, so the appropriate velocity functional may be also normalized by the 
total length. 

To simplify the analysis we further assume that the motion of the active segment is steady [32, 33] with unknown 
velocity V = L = /. and that the unknown functions <x, v and the unknown controls £,x depend exclusively on the 
appropriately chosen co-moving coordinate u - (x - xq - Vt)/Lo e [-1/2,1/2] . Then in dimensionless variables 
cr := cr/x* , x := xj yjtj/i;*, t := t/(tj/x*), % '■= £/£* and;^ '■= x/X* we obtain the force balance equation 


1 a / d u cr(u) 

i 2 "l *t(«) 


• oi u) = g 2 (u). 


( 6 ) 


where L Lq/ y/rj/t;*. The re-scaled control functions 

gi(u) = £(Lu) > 0,g 2 (u) =x ( Lu ) ^ 0 


must satisfy the constraints 


The boundary conditions take the form 


where 


nU2 nl/2 

I gi(u)du - I g 2 (u)du = 1. 

J-l/2 J-i 


- 1/2 


cr(-l/2) = cr(l/2) 

j_ ( d„o\ 1/2) _ a„cr(-i/2) \ _ _-rr7 

L2 ( ifid/2) gi(-l/2) / - 


(7) 


( 8 ) 


AV:= (V+ ~ V - )n 


The dimensionless velocity of the segment per length V = V/L can be found from the formula 

T7_ — , 1 / d»°~(l/2) , d u( r(-H2) \ 

2L 2 \ gl ( 1/2) gi(-l/2) ) 


( 9 ) 


— _ (v+ + V-)tj 

* m • — 


where 


2t*T 0 
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Now, if we assume that the two parameters ( V m ,AV ), characterizing actin treadmilling, are fixed we can pose the 
optimization problem of finding the controls gi(u), giiu) ensuring the maximization of the normalized velocity V. 
This problem is nontrivial because the functional V{g\,g 2 ] is prescribed implicitly through the unknown solution 
of the boundary value problem ( 6 , 8 ). To our advantage though this linear elliptic problem is classical and is well 
understood, e.g. [54], 

We observe that parameter V m enters the expression for the velocity (9) in an additive way and does not affect 
the solution of the optimization problem. The reason is that V m characterizes a propulsion mode associated with 
simple accretion of the material at the front and its simultaneous removal at the rear; when V m 40ana priori polarity 
is imposed and the problem of motility initiation disappears. In view of the complete decoupling of this mode of 
self-propulsion from our controls, in what follows we assume without loss of generality that V m = 0. 

The second parameter AV, also characterizing the protrusion strength, does not induce polarity. As we clarify 
in the next Section, this parameter represents the mechanical action of pushers and the dependence of the crawling 
velocity on AV, which is now sensitive to both controls, is much more subtle than in the case of V m . 

3. Pushers and pullers 

To see that parameter AV characterizes in our setting the activity of pushers, consider the global balance of couples 
in the co-moving coordinate system 

/-1/2 _ n 1/2 

L I g\(u)v(u)udu — -AV + I g 2 (u)du + (Tq. (10) 

J-l/2 J-l/2 

Here the term in the left hand side characterizes the total moment due to external (frictional) forces. The first term in 
the the right hand side 

T = AV 

f 1/2 

is due to active protrusion, while the second term J_ | 0 g 2 (u)du - 1 is due to active contraction. The last term cr 0 
corresponds to passive reaction forces resulting from the prescription of the length of the segment. 

Our assumption that AV > 0 means that the protrusion couple has a negative sign showing that the corresponding 
force dipoles act on the surrounding medium by pushing outward and creating negative stress. Instead, the contraction 
couple has a positive sign because the contractile forces pull inward and the induced stresses are positive. We can 
therefore associate protrusion with the presence of pushers and contraction with the activity of pullers [18, 21]. We 
can also (tentatively) argue that motility is protrusion-dominated when T > 1 and it is contraction-dominated when 
0 < T <1. This assertion will be confirmed later in the paper by rigorous analysis. 

To illustrate the different roles played in our motility mechanism by pushers and pullers, we present below an 
analysis of a toy model which anticipates the main conclusions of the paper. We temporarily set AV = 0 and describe 
the distribution of pushers and pullers by the same function gi{u) allowing it now to be both positive and negative. 
Our goal is to show that for protrusion-dominated motility driven by pushers, it is beneficial to create strong adhesion 
at the leading edge while for contraction-dominated motility driven by pullers, it is the trailing edge that has to adhere 
most strongly. 

Consider a special choices of control functions, 

g\(u) = q6(u - U\) + (1 - q)6(u - i<2),g2(u) = pd(u - M3) - (1 - p)6(u - M4), (11) 

where 6 is the Dirac distribution, 0 < q,P< 1 and -1/2 < u\, U 2 , M 3 , M 4 < 1 /2. The control function g\(u) represents 
two adhesion sites u — u\ and u = U 2 whose locations and intensities must be optimized in relation to the prescribed 
position of a single puller placed at n = M 3 and characterized by a positive dipole moment p and a single pusher located 
at m — M4 and characterized by a negative dipole moment p — 1 . The parameter p measures the relative importance of 
contraction comparing to protrusion. 

Suppose that 

M3 < Ml, M2 < M4 (12) 
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which ensures that our active segment moves from left to right and that in the segment the adhesion sites are always 
outside the location of the active agents. By using (11) we can express the velocity of the segment as 

y _ 1 L + (1 - 2 p)(uz - Ui)c \ 

2 \ u\ - ui + 1 + (ui - ii\)a) ’ 

where, 

1 4 1 - 2 q 

a = -- and c = --. 

1 + q(\ - q)(ui - U\)L - 1 + q(\ - q)(ui - U\)L~ 

Suppose for simplicity that our two adhesive complexes are placed symmetrically with respect to the center of the 
segment u\ = 1 — no. We can then use a single parameter Am — u\ + 1/2 = 1/2 - ui to obtain, 

V = I (i + (1 - 2 p)(l - 2 q) 

2 \+ 2 Auq(\-q)L^ + ^ u / 

At Au = 0, when adhesive complexes localize at the edges of the segment, the velocity reaches its maximum value 

V= l -[\+(\- 2 p)(\- 2 q)\. (13) 

Because of the imposed inequalities (12), we must necessarily have in this configuration M 3 = -1/2 and M 4 = 1/2 
meaning that pushers and pullers must also localize at the edges. Notice however that inequalities M3 < u\ and ui < 114 
are necessary for formula (13) to hold so the apparent conclusion that M3 = u\ — 1/2 and 114 = 112 = -1/2 is a result 
of an abuse of notation. In a theory where, pushers, pullers and adhesion complexes have a characteristic size of 
dispersion, the adhesion clusters will be slightly ahead of pullers and slightly behind the pushers so that active agents 
take advantage of the firm attachment to either push or pull (see more on this subject below). 

From our simple analysis it follows that if pullers dominate (p > 1 /2), the optimal way to tune adhesion is to set 
< 7=1 and concentrate adhesive complexes at the trailing edge of the moving segment achieving maximum velocity 
V = 1 when p = 1. If, instead, pushers dominate (p < 1/2), the optimal way to tune adhesion is to set q = 0 and 
concentrate adhesive complexes at the leading edge. In this case, the maximal velocity is again V = 1 given that 
p = 0. In other words, to be effective, pullers have to concentrate on the trailing edge and ensure strong adhesion on 
the leading edge: in this way pullers can inflict contraction that displaces the trailing edge which due to the length 
constraint also propels the leading edge. On the contrary, pushers can take advantage of the firm attachment at the 
trailing edge to push against it and propel the leading edge which in turn pulls the trailing edge due to the length 
constraint. 

4. Contraction driven motility 

We now return to the study of the optimization problem in the original formulation. The simplest 
transparent case is when protrusion is disabled AV = 0 and motility is fully contraction-driven. 

Suppose first that g\ = 1 which means that the adhesion complexes are distributed uniformly. Then 
of the active segment can be expressed as a quadrature 

1 r 1/2 

V = ~ ■ , 77 r sinh (Ln)g 2 (u)du. ( 14 ) 

2 stnh(|) J- 1/2 

One can see that if the function gi(u) is even, then V = 0. This result can be interpreted as an analogue of the famous 
Purcell’s theorem about the necessity of non-symmetric strokes in Stokes swimming [4, 2]. If the distribution gi(u) is 
non-symmetric and, for instance, more motors are placed at the rear of the segment, the velocity will become positive. 
Using the fact gi(u) > Owe can also conclude from (14) that V < 1/2. This upper bound is reached when all the 
motors are fully localized at the rear and gi(u) = 6(11 + 1 / 2 ). 


analytically 
the velocity 
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Now, consider the general case when the focal adhesions are distributed inhomogeneously. Since (6) is a Sturm- 
Liouville problem, its solution can be written as 


ct(m) = cro - 



G(m, i)[g 2 (i)-cr 0 ]c/i. 


(15) 


where the Green’s function G(u, s) can be represented two auxiliary functions hiu) and f{u) 

G(u, s) = ^ [h(u)f(sW [s<u] + h{s)f{u )¥ [u<s] ], 

solving the following standard boundary value problems [54] : 

f {—h'y = L 2 h ( (i-/')'=L 2 / 

{ h(-1 /2) = l,/t(l/2) = 1 ’ { /(-1/2)‘= l,/(l/2) = -1 ' 


(16) 


(17) 


In (16), C = (hf - fh')lg\ is a constant involving the Wronskian of the two auxiliary functions h(u) and f(u) and Ik 
is the indicator function. We can now write 

- 1 f 1/2 

V =T f(u)(giiu) - g 2 )du, (18) 

1 J-1/2 

where we introduced a new measure of inhomogeneity of contraction: 

XV/2 h(u)g2iu)du 

g -■ — .. | j 2 • 

J_ l j2 hiu)du 

If both functions g i, 2 (m) are even, then fin) is odd and, since the integral of a product of an odd and an even functions 
is equal to zero, we obtain that V — 0. The same result follows if we assume that contraction is homogeneous 
g 2 (u) = £2 = 1 while the adhesion distribution g\iu) is arbitrary. Therefore, to ensure motility at AV = 0, contraction 
must be inhomogeneous while adhesion may still be uniform provided contraction is not even. 

To find the optimal distributions giiu), giiu) we proceed in two steps. We first show that V < 1 and then find a 
configuration of controls allowing the cell to reach this bound. 

Notice that we can rewrite (18) in the form 


v = \ ( J' f{u)ig 2 (u) 



f{u){g 2 {u) - gi)du 


where we defined the domains S _ = { u/g 2 iu ) < g 2 ] and S + = {u/g 2 iu) > • Applying the maximum principle to 

(17) we obtain that 1 > hiu) > 0 and hiu ) > fin) > -hiu). Using the bounds on /, we can write 


V < 


m 


h(u)g 2 iu)du + g 2 I hiu)du 


L 


Since the integrands are positive and liiu) < 1 it finally follows that 


_ pl/2 rW 

V < I hiu)g 2 (u)du < I g 2 iu)du = 1. 
J-1/2 J-1/2 


(19) 


Observe that in the case of a homogeneous distribution of adhesive clusters, the velocity could reach only one half of 
this maximal value. 

We now show that the maximal velocity V = 1 can be reached if both controls g\(u) and g 2 iu) are fully localized. 
Take 6 > 0 and consider a regularized distribution 


giiu-,6) 


1 0 
n 6 1 + iu - «i) 2 
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For this choice of gi(u) the auxiliary functions h(u) and f(u ) can be written explicitly in term of Legendre polynomials. 
In the limit 0 —» 0 and lim^o giU< ; 0) = 5{u - ni) we obtain 


h(u) = 1 and f(u) = { ' ^ U Ul 

J [ -1 if u > U\. 

By using these explicit expressions for the auxiliary functions we can rewrite (18) in the form 


if we now suppose that g 2 (u ) = 6 (u - u 2 ) the expression for velocity reduces to reduces to 

j ( 1 - 2 u\ if «2 < Mi 

V = - < —2mi if M 2 = Mi 

( — 1 — 2 Mi if 142 > Ml 

It is now clear that the velocity reaches its maximal value as Mi —» -1/2 while M 2 < mi. We can then formally write 
U 2 = Mi = -1/2 and obtain the controls g 2 (u) = gi(u) = 6 {u + 1/2) saturating the bound V = 1. Notice, however, that 
if we assume directly u\ = M 2 —» -1/2 in (20), we obtain V — 1/2. This is in agreement with the physical intuition 
that the anchorage point must be located to the right of the pulling force dipole: in this case the pulling forces advance 
the rear edge of the segment with minimal slipping. Mathematically, we encounter here the case of non-commutation 
of the limiting procedures M 2 —> -1/2 and U\ —> -1/2 giving V = 1 only if the limits are taken in the above order. 

To summarize, the optimization of the distribution of focal adhesions allows the contraction-dominated mecha¬ 
nism of cell motility to reach the value of velocity which is twice as large as when the adhesion is uniform. This 
means that in order to improve the motility performance the adhesion machinery must conspire with the contraction 
machinery making sure that both the motors and the adhesive centers are localized at the trailing edge. Interestingly, 
exactly this type of correlation between the stresses created by contraction and the distribution of focal adhesions was 
observed in experiments and numerical simulations [55, 56, 57, 58, 59, 22], The localization of adhesion complexes 
close to cell edges, where contraction is the strongest, has been also reported outside the motility context [60, 61, 62], 


X M1 /'■!/•£ 

g2(u)du - I g2(u)du - 2 mi 

1/2 Jh, 


( 20 ) 


5. Upper and lower bounds for velocity in the general case 


We now turn to the general case where both contraction and protrusion are active. In particular, the protrusive 
power will be characterized by the parameter AT - T > 0 which was assumed to be equal to zero in the previous 
Section. We can then write 



ri'mdu r m 

- t + m( 82 (u)-g2)du 

f 1/2 h(u)du J- 1/2 


( 21 ) 


As we see, the first term in the right hand side is associated with protrusion-based (or filament-driven [63]) motility 
while, as we have already seen, the second term is the contribution due to contraction-based (or motor-driven [63] ) 
motility. We notice that if gi(u) is even, then /(n) is odd and h(u) is even, leading to 


XV/2 /(«)^M 

j ^ h(u)du 


If g 2 (u) is also even, then 



f(u)(g 2 (i4) - g 2 )du = 0. 


In this case the velocity of the segment is fully controlled by the accretion mechanism characterized by the parameter 

VV 
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Consider first the case of protrusion-driven motility by assuming that contraction is homogeneous g 2 (u) = 1 and 
therefore does not contribute to the overall velocity. By using again the maximum principle we obtain inequalities 


J.'i^ f( u )du 

XV/2 h ( u ) du 


< 1 , 


leading to the upper bound 

V=Vp< T -. (22) 

The maximum of the protrusive contribution to velocity is reached when, g\(u) = d(it - y), because in this case h = 1 
and / = 1 almost everywhere. Observe, that the optimal solution in the case of protrusion-driven motility is in some 
sense opposite to the solution gi(u) = 6{u + 1 / 2 ) obtained in the case of the contraction-driven motility. 

Based on (19) and (22) we can now argue that in the case when both treadmilling and contraction are present, an 
upper bound for velocity is 

- T 
V < - + 1, 

2 

however, as we have just seen, in view of the incompatibility of the corresponding optimal controls, this bound cannot 
be reached. The optimal strategy for focal adhesions would then involve a compromise between the necessity to 
localize adhesion at the trailing edge in order to assist the contraction mechanism and the competing trend to localize 
adhesion at the leading edge in order to improve the protrusion power of the cell. 

To obtain a lower bound for V we now consider a particular test function representing a weighted sum of our 
competing optimal controls, gi(u) = qd(u + 1/2) + (1 - q)6(u - 1/2). We also chose g 2 {u) = S{u - mq), where q 6 [0,1] 
and uq e [-1/2,1/2] are two parameters to be optimized. Then, by solving (17) we obtain. 


and, 


f(u) = 


h(u ) 


1 if u = -1/2 

2 if u e] - I/! 
-1 ifn = 1/2 


TTgjFif" e]-l/2,l/2[ 


1 ifw = -1/2 
p- if m e] - 1/ 
1 if u = 1/2, 


TTin L_if«6]-l/2,l/2[ 


which leads to the expression for the velocity 


V = ^-(1 - 2 q) + i(/(« 0 ) - (1 - 2 q)h(uo)). 


The optimization with respect to m (i gives m (i = -1/2 and 


V=--q(T- 1). 


Finally, optimizing in q we obtain that if T < 1, we need to take q - 0 and if T > 1, we get q = 1. This result, illustrated 
in Fig 1, suggests that there is a switch at T — 1 between the contraction-centered optimization strategy (q = 0) and 
the protrusion-centered optimization strategy (q = 1). Notice that the switch takes place exactly when the negative 
protrusion generated couple T becomes equal to the positive contractile couple equal to 1. At an interesting state T = 
1 , the two active mechanisms neutralize each other and activity-related dipoles become invisible behind the passive 
contributions in Eq. (10): in this case the optimal position of active and adhesive agents becomes indeterminate. 
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Figure 1: Solid lines: Lower bound on the optimal velocity of self-propulsion V as a function of the measure of the (relative) protrusive strength 
T. The optimal strategy depends on whether contraction (T < 1 ) or protrusion (T > 1) dominates. The dashed line represents the upper bound 
obtained by formally summing the incompatible upper bounds for the protrusion and contraction based strategies. The dotted line represents a 
sub-optimal strategy obtained under the assumption that adhesion is homogeneous. Insets illustrate the associated configurations of controls g\{u ) 
and g 2 (u). 


6. Numerical solution of the optimization problem 

To show that the low bound obtained in the previous Section is rather close to being optimal, here we solve the 
optimization problem numerically. A finite dimensional reduction of the original variational problem is constructed 
by selecting N + 2 points u, = i/(N + 1) - 1/2 that subdivide the segment [-1/2,1/2], We then localize adhesion and 
contraction in these points by choosing the control functions in the form 

N N +1 

8 1 («) = ^g'lSiu - Ui),g 2 (u ) = ^g‘ 2 6(u - ut), 

1=1 1=0 


where E/Io* 8 o ~ 2/Ii 8 \ = 1 and g\ > 0, g ' 2 > 0. In this way we also mimic the discrete nature of real adhesive 
clusters and real myosin motors, so in some respects the discrete problem is even more realistic than the continuum 
one. 

By solving the auxiliary linear elliptic problems for this choice of controls we obtain that the functions h(u) and 
/(m) are piece-wise constant h(u) = [«;.«, an d /(«) = E/Io C^[k ( T he coefficients with i e [2, A] 

satisfy the equations 

s'lCVt - A,_ 2 ) + g\gX l Ai^L 2 { Ui - Uj -,) = g‘[ l (Ai - A,-_i) 

£i(C,-t - Ct- 2 ) + g[g'^ l Ci-iL 2 (ui - «,--!) = g\-\Ci - Ci_i) 

The boundary conditions give 

Aq = 1,A^ = l,Co = 1 ,C N = -1. 


We use the conventions A [ = Aq, A N+ \ = An, C \ = Co and C ; v+1 = C,v to express the velocity V in the form 



N +1 

r-Z 

i =0 


i Ai + A,_] 
8 2 ^ 


ZZo C i( U i +1 ~ “/) 


1 V ,■ C«- + Cj-1 
+ 2Zj ?! 2 

i=0 


(23) 
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The function (23) was optimized numerically with respect to parameters g l { and g‘ 0 subjected to the appropriate equality 
and inequality type constraints. To find the global minimum we used the method of simulated annealing with initial 
guesses corresponding to homogeneous configurations. In Fig. 2, illustrating our results for N = 100, one can see that 
for T < 1 both optimal functions g\(u) and giiu) are localized at the trailing edge. Instead, for T > 1 we observe 
that gi(u) localizes at the leading edge while g 2 (w) localizes at the trailing edge. The value of the maximal velocity, 
obtained numerically, is the same as in the bound up to an error proportional to the mesh size. Notice also that the 
optimally spaced adhesion points are shifted from the locations of force dipoles by a mesh size. Thus, the solution 
of the numerical (regularized) problem in the contraction-dominated regime (with T = 0.9), shown in Fig. 2, clearly 
distinguishes the optimal functions g\{u) and y 2 00 that are both localized at the size of the mesh. The function 
g\{u) remains different from zero at one mesh size beyond the point where we already have giiu) = 0 (for positive 
velocity). In the protrusion-dominated regime (at T = 1.1) shown in Fig. 2, the mesh size again prevents localization 
of the function g\{u) exactly at the leading edge while g 2 («) localizes exactly at the trailing edge. These observations 
confirm that the regularized numerical solution is more realistic than its singular analytic prototype. Overall, our 
numerical results are in full agreement with the analytic bounds which suggests that that these bounds are (nearly) 
sharp. 



Figure 2: Numerical results for the optimization of the function (23) for two values of the protrusive strength, below the threshold T = 0.9 and 
above the threshold T = 1.1. Insets magnify the structure of the optimal controls near the boundaries of the moving segment. Parameter L = 10. 


7. Local stability analysis 


In this Section we use a perturbation analysis to provide additional evidence that our lower bounds are close to 
being optimal. 

In what follows we use the superscript o to indicate the unperturbed state and superscript ★ to mark parameters 
characterizing the perturbed configuration. We assume that for all functions <p{u) the following expansion holds in the 
first order 

<p(u) — cf>°(u ) + e<p*(u ), 

where e is a small parameter. Keeping only the first order term in the expression for V, we obtain 


-★ 

V 


where h* = H*\ f* = F*' and 


l 


m sT 


+ f°sl 

T ~f-W2 h °«°21 r 1 / 2 

fl /2 I a 

J-W 2 h 


f 1 ' 2 r 

J—1/2 J 

r 


(f- 




I fl/2 _ Q* f 1/2 fo \ 

\ J - 1/2/ fZ h ° J - 1/2f ) 



H*" - L~g\H * = L~h°g* ( F*" - I?g\F * = I?F°g * 
H*'(- 1/2) = H*'( 1/2) = 0 ’ \ F*'(- 1/2) = F*'(l/2) = 0. 


( 24 ) 


In view of the saturation of the constraints by the unperturbed solution 


rZ/2 S2 = fV/2 we mUSt have 



r 1/2 

sf = 0 . 

J-l/2 
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Suppose now that we perturb the ’optimal’ controls g°(u ) = g°{u) = 6(u + 1/2) delivering the lower bound for 
velocity at T <1. We can again formally set 


h°(u) = 1 and/°(i<) = j 


1 if m = - 1/2 
if u > - 1 / 2 , 


while remembering that the localization point for adhesion must be shifted with respect to the point of localization of 
contraction. The perturbation of velocity can be written as 


-* 1 
V =2 


r l/2 pi/2 

(/» + h°(i,))g*(u)d u + (T - 1) (/*(«) + h 

J-l/2 J-l/2 


(u))du 


A rather general class of perturbed controls can be represented in the form 

j g*(u ) = -q5{u + 5 ) + r(u) with, r(u)du = q 

| g*(u) = -pS(u + j) + l(u) with, j' ] ^ l(u)du = p, 

Since gi > 0 and g 2 > Owe demand that r(u) > 0 and l(u ) > 0 and therefore also q > 0 and p > 0. We obtain 


V* = -p + (T-l)£ 


1/2 


r(u)(u + 1 /2)du. 


(25) 


(26) 


1/2 


From (26) we see that if T < 1, a perturbation of the controls g^u) = g° 2 (u) — 5{u + 1/2) leads to the decrease of the 
-★ 

velocity: V <0. Instead, if T > 1, by choosing r(n) such that. 


rm 

1 

J-l/2 


1 p 

r(u)(u H —)du > 


t - r 


we obtain a perturbation with V >0. These observations suggest that beyond the threshold T — 1 the control 
function g\(u) should be no longer localized at the trailing edge edge. We also see that the most ’efficient’ way to 
make the velocity larger at T > 1 is to localize the function r(u) at the leading edge of the segment (at it - 1 /2). 

Consider now a perturbation of the set of controls g\(u) = 6{u - 1 /2) and g°(u) = 6{u + 1 /2) delivering our lower 
bound for T > 1. With the same abuse of notations as before we can write 


h°(u) = 1 and f°(u) = 


1 if u < 1/2 
-1 if u = 1 / 2 . 


We represent the perturbations in the form 


{ g*(u) = —qd(u ~ {) + r(u) with, r(u)du = q > 0 
g*(j/) = -p6(u + 5 ) + l(u ) with, l(u)du = p > 0 . 


(27) 


1/2 

where again r(u ) > 0 and I(u) > 0 and therefore also q > 0 and p > 0. The ensuing perturbation of velocity is 


-* 1 
V =2 


r l/2 

( 

J-l/2 


(f°(u) - h°(u))gl(u)du + (T- 1) (f*(u) - h (u))du 


r i/2 

( 

J-l/2 


= -(T- 1) 


X 


1/2 j 

i/2 ( 2 


u)r(u)du 


It is now clear that if T > 1, then V* < 0, showing that the perturbations of controls are sub-optimal. This gives 
another evidence that the test function providing the lower bound for velocity at T > 1 are at least close to being 
optimal. 

Based on this analysis we conjecture that the function V(T), representing the optimal velocity, is piece-wise linear 
with a kink at T = 1. The presence of a threshold indicates a switch from contraction-dominated motility pattern 
to protrusion-dominated motility pattern. As the relative power of protrusion, epitomized by T, increases beyond 
this threshold, the focal adhesions, maintaining the optimality of the self-propulsion velocity, must migrate from the 
trailing to the leading edge of the active segment. The dynamic migration of adhesion proteins to the edges has been 
observed in experiments [60]. In real cells, however, both edges are usually populated by adhesion complexes and we 
can speculate that in this way cells can adjust smoothly to transitions from one driving mode to another. 
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8. Discussion 


In this paper we used a simple analytically tractable model of cell motility to study the optimal strategies allowing 
cells to move faster by actively coordinating spatial distributions of contractile and adhesive agents. Our study reveals 
that if adhesion complexes can detect the dominating mechanism of self propulsion, they can self-organize to ensure 
the best performance. 

We made specific predictions regarding the advantageous correlations between the distributions of adhesive and 
force producing agents and showed that the dependence of the maximal velocity of self-propulsion on the relative 
strength of contraction and protrusion may be non-monotone. In particular, our model predicts that a limited activation 
of protrusion will necessarily lower the maximal velocity achieved in a purely contractile mode of self-propulsion. 
However, as the protrusion strength increases, protrusion can overtake contraction and the velocity of self-propulsion 
will increase beyond the level achieved in the contraction-dominated case. 

Previously we have shown that contraction-driven motility mechanism may be sufficient by itself to explain cell 
polarization, motility initiation, motility arrest and the symmetrization of a cell before mitosis [31, 44]. From the 
analysis presented in this paper it becomes evident that, if the speed of self propulsion is an issue, cells should 
mostly rely on protrusion. More specifically, to maximize its velocity performance after motility initiation a cell must 
switch from contraction-dominated to protrusion-dominated motility mechanism by increasing the protrusive power 
and appropriately rearranging the distribution of adhesive complexes. It was shown in [43] that similar transitions 
between contraction and protrusion mechanisms can be used by a cell to accommodate different types of cargos. 

To compare our predictions with experiments we can use numerical values of parameters for keratocyte fragments 
[29, 34, 32]. We obtain the following rough estimates: x* = 10 3 Pa, f = 3 x 10 16 Pfl ■ m -2 ■ s, r\ = 3 x 10 4 Pa • s, 
Lo = 10 inn, V m = tiftm ■ min 1 and AV = 0.6 inn ■ min -1 . Our first quantitative prediction concerns the case when 
active treadmilling of actin is knocked down and adhesion is homogeneous. As we have shown, the largest velocity 
in this case is reached when all myosin motors are localized at the trailing edge as observed in most eukariotic cells 
[64, 65, 66, 67]. In dimensional form, the predicted maximal velocity is V — L^x"'/(2rj) « 10 fim ■ min -1 , which is low 
in view of the data on keratocyte fragments suggesting that velocity should be in the interval 30 - 40 fim ■ min -1 [68]. 
This is not surprising because according to our results only half of the total amount of motors is ’’used” in this case 
which implies that adhesion homogeneity is highly sub-optimal. 

If the adhesion inhomogeneity is allowed, the configuration becomes optimal when both myosins and integrins are 
localized at the trailing edge. Such highly correlated distributions have been observed in experiments and generated 
in microscale-based numerical models [60, 56, 61, 62, 57, 58, 59, 69]. In this case, as a result of the cooperative 
response, a cell can be more efficient achieving the velocity that is two times larger than in the case of homogeneous 
adhesion: V = Lqx* h ~ 20 inn ■ min -1 . 

To be even closer to reality we need to take active treadmilling into consideration and our estimates suggest that 
T ~ 0.1 <sc 1. This means that we are in the contraction-dominated motility regime for which our velocity bound 
gives more realistic value V = V m + Lqx* h ~ 28 jim ■ min -1 . Notice, however, that reducing the valuer* by one order 
of magnitude, which is within the existing error bounds, we may easily reach the regime where T > 1 and where 
it would becomes more beneficial for adhesion clusters to localize at the leading edge conspiring with protrusive 
elements. A spatial correlation of this type between adhesion and protrusion has been recorded in both experiments 
and comprehensive numerical models, see for instance the data in [70, 71] on nematode spermatozoa. 

The proposed model can be also tested indirectly. For instance, we know that the location of adhesion complexes 
in a moving cell can be identified by measuring the distribution of traction forces in the elastic environment [72, 73]. 
If the adhesive complexes are found to be shifted towards the leading edge, we would argue that the cell relies for its 
advance mostly on actin treadmilling. If instead the adhesive complexes are preferentially positioned at the trailing 
edge, our model suggests that motility is mostly driven by contraction. Both predictions can be tested by independent 
measurements. 

An interesting possibility would be if cells could alternate the location of maximum adhesion between the trailing 
edge and the leading edge in response to oscillations in the level of activity of pullers and pushers. The evidence of 
such switching may be that both location are populated with adhesive complexes. We may also recall that the classical 
mechanism of crawling for eukaryotic cells involves two phases [24, 74, 26, 75]: one of them is associated with the 
creation of protrusions that push at the front while relying on the stabilization of the trailing edge and another one 
involves pulling of the rear (of the cargo) which requires fixation at the front. The switching between these two phases 
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takes place almost periodically and the associated reorganization of adhesion clusters from the leading to the trailing 
edge has been well documented [76, 72]. To capture such non-steady motility pattern in the framework of our model, 
the simplifying traveling wave assumption would have to be replaced by a more complex ansatz. 

In conclusion, we emphasize that our interpretations are based on the study of a variational problem whose analysis 
revealed some interesting correlations between the spatial arrangement of adhesion and contraction agents and has 
led us to quantitative predictions that are in agreement with experiment. The prototypical nature of the proposed 
model, however, conceals considerable complexity of the actual cell motility phenomenon which involves intricate 
bio-chemical feedback loops, geometrically complex mechanical flows and highly nontrivial rheological behavior. In 
particular, the singular nature of the obtained solutions can be at least partially linked to the fact that treadmilling 
is modeled schematically, with polymerization and depolymerization processes localized at the edges: at least one 
additional control function describing the distribution of pushers is needed to regularize the problem in this respect. 
The situation is complicated further by the fact that the dominant trade-off condition, controlling the self-organization 
of active agents, is still unknown notwithstanding some recent results in this direction [23], However, even in the 
absence of the definitive optimization criterion and with minimal assumptions about the inner working of the motility 
machinery, our study reveals that depending on the task and the available resources a cell may have to modify its mode 
of operation rather drastically to ensure the best possible performance. 
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